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SFCTION  I 


INTRODUCTION 

For  the  case  of  flow  over  an  airfoil  in  a freestream  at 
Mach  numbers  near  1,  small  amplitude  motions  of  the  body  surface 
can  produce  larqe  variations  in  the  aerodynamic  forces  and  moments 
actina  on  the  structure.  In  addition,  phase  differences  between 
the  flow  variables  and  the  resultant  forces  may  be  areat.  These 
characteristics  tend  to  enhance  the  probability  of  encounterinq 
aeroelastic  instabilities  in  the  transonic  flow  reqime  and  thus 
evidence  a need  for  techniques  of  analyzing  both  the  flow  field 
and  the  structural  response  for  such  situations. 

In  the  subsonic  or  supersonic  case,  the  leadina  order  flow 
equations  are  linear  such  that  the  aerodynamic  forces  depend 
upon  the  body  motion  in  a linear  fashion.  Furthermore,  the 
resultant  forces  actinq  on  the  airfoil  may  be  obtained  throuqh 
superposition  by  summing  the  contr ibut ions  due  to  the  various 
types  of  body  motion  being  considered.  This  allows  the  linear 
structural  equations  of  motion  to  be  solved  independently  of  the 
aerodynamic  equations  which  provide  only  the  force  coefficients. 
Uncouplinq  of  the  fluid  and  structural  equations  is  not,  in 
general,  possible  for  the  transonic  reqime  due  to  its  inherent 
nonlinear  nature. 

Advances  in  computational  methods  (1-16)  have  made  a number 
of  techniques  available  for  computing  unsteady  transonic  flows. 
While  several  different  physical  problems  have  been  considered, 
the  unsteady  body  motion  was  generally  prescribed  as  a known 
function  of  time  thereby  precludinq  the  simulation  of  aeroelastic 
behavior.  It  is  only  recently  that  the  application  of  these 
computational  procedures  to  actual  aeroelastic  problems  has 
appeared  (17,18). 

It  is  the  purpose  of  this  report  to  describe  a method  for 
analyzinq  the  structural  response  of  a two-dimensional  airfoil 
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in  transonic  flow,  and  to  provide  computational  examples  of  the 
results  obtained  by  applying  this  procedure  to  physical  situations 
of  practical  interest.  The  governing  aerodynamic  equation  of 
motion  is  assumed  to  be  the  unsteady  low  frequency  small 
disturbance  transonic  equation  for  the  velocity  potential  func- 
tion which  is  capable  of  simulating  nonlinear  flow  phenomena 
including  irregular  shock  wave  motions.  Structural  equations 
of  motion  are  formulated  for  a three-degree  of  freedom  airfoil 
by  modelling  the  structure  as  a spring-mass  system.  The  coupled 
aerodynamic-structural  equations  are  then  simultaneously  inteqrated 
in  time  such  that  the  flow  field  and  the  response  of  the  airfoil 
to  tne  resultant  aerodynamic  forces  are  allowed  to  interact  in 
a manner  much  like  the  physical  situation. 

The  method  of  time  integration  has  already  been  applied  in 
a superficial  manner  for  analyzing  a one  decree  of  freedom 
airfoil  (17).  Tn  this  previous  work,  several  types  of  motion 
were  produced  by  varyinq  the  structural  parameters.  An  alternative 
point  of  view  is  taken  in  this  report.  Here,  the  structural 
parameters  are  presumed  fixed  and  the  motion  resulting  from 
various  choices  of  the  initial  conditions  is  considered.  The 
airfoil  selected  for  study  is  an  NACA  64A010  airfoil  which  is 
10?  thick  and  representative  of  transonic  airfoils  currently  in 
use.  Plunginq,  airfoil  pitchinq,  and  aileron  rotation  decrees 
of  freedom  have  been  allowed  and  structural  parameters  have  been 
assigned  representative  values.  Solutions  were  obtained  for 
several  choices  of  initial  conditions  and  non-dimensional 
freostream  density.  The  effects  of  these  parameters  on  the 
time-dependent  structural  response  and  aerodynamic  coefficients 
is  presented. 
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The  system  (23)  - (2S)  may  be  written  as  the  vector  equation 
MX"  (t)  + PX'  (t)  + KX(t)  = r (r)/u 


(26) 


where  the  initial  conditions  X(0)  and  X'(0)  complete  a description 
of  the  problem. 

2.2  Unsteady  Transonic  Potential  Equation 

A number  of  procedures  are  available  for  the  solution  of 
unsteady  transonic  flow  problems.  Methods  of  harmonic  analysis  (3-7) 
assume  that  the  unsteady  flow  may  be  considered  as  a small  perturba- 
tion about  a steady  state  due  to  the  motion  of  the  boundary.  These 
time-linearized  methods  fail  to  treat  shock  motions  exactly  and  thus 
are  limited  in  application  to  very  small  amplitude  airfoil  motions 
which  constrain  shock  waves  to  small  displacement  about  their  steady- 
state  position.  In  addition,  the  unsteady  perturbations  are  prescribed 
to  have  a harmonic  dependence  in  time.  While  such  methods  are  not 
capable  of  treatinq  qeneral  airfoil  motions,  they  can  provide  input 
for  classical  aeroelastic  computations  (18). 

A more  exact  treatment  of  unsteady  transonic  flow  problems  is 
provided  by  integratinq  the  equations  of  motion  in  time.  Solutions 
to  the  unsteady  Euler  equations  (8,8)  and  the  full  potent  ial 
equation  (10)  have  been  obtained  by  explicit  finite  difference 
schemes.  Althouqh  these  methods  can  in  principle  be  applied  to 
quite  general  airfoil  motions,  they  are  limited  bv  a time  step 
restriction  for  eomputat lonal  stability  thereby  making  practical 
calculations  prohibitive.  If  consideration  is  limited  to  irrota- 
tional  flows  and  low  frequency  motions,  then  a reduced  form  of  the 
unsteady  potential  equation  results.  A fully  time-implicit 
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method  of  solution  of  this  equation  has  been  developed  (17)  such 
that  no  time  step  restriction  for  stability  exists.  In  addition, 
this  equation  is  capable  of  simulatinq  nonlinear  flow  phenomena 
includinq  irregular  shock  wave  motions.  Solutions  of  the  unsteadv 
low  frequency  potential  equation  have  compared  nuite  well  with 
solutions  of  the  Euler  equations  (17) . 

We  now  consider  the  two-dimensional  irrotutional  unsteady 
flow  over  an  airfoil.  The  velocity  potential  function,!  , is 
expanded  as 

( x , y , t ) = umc62 /3<J>  (f.  ,n,T)  (27) 

where 

n = yfl  ^ 3 /c , (29,) 

Z and  t are  given  by  equations  (9)  and  (10)  respectively,  and  it 
is  assumed  that  6<<1.  If  these  expressions  are  substituted  into 
the  full  unsteady  potential  equation  and  the  coefficients  of  like 
powers  of  are  equated,  then  the  leadino  order  result  is 


[(1  - M2 ) 6~ 2 / 3 

U CO 


(1  + y ) M 2 

1 <x>  r r r 


+ * 


(28) 


if  it  is  assumed  that  (1  - M2)^-2/3  = 0(1). 

This  is  the  form  of  the  unsteady  low  freouenev  small  disturbance 
transonic  potential  equation  which  will  be  considered  appropriate 
for  describinci  the  flow.  Correspondinq  boundary  conditions  for 
equation  (29)  are  obtained  from  the  flow  tanqencv  condition  on 
the  airfoil  surface,  from  the  Kutta  condition  at  the  trail  inn  edoe 
with  a constant  jump  in  potential  across  the  vortex  sheet  in  the 
wake,  and  by  requiring  the  perturbation  velocity  to  vanish  far  from 
the  body. 

The  unsteady  airfoil  boundary  is  defined  by 


y - F(x,t)  = 0 for  0 £ x < c. 


( 10) 
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Hero 

F(x,t)  [ f ( ) + o(v)A"!  - (;.  - •o)>'»(t) /«$  - (•  - f )!l(f.  - .:.f)t<(i)  S]  HI) 

whore  f ( ;.)  corresponds  to  the  specified  airfoil  geometry:  a , a , and 
n are  the  previously  defined  plunging,  airfoil  pit  china,  and  aileron 
pitching  displacements  respectively;  and  !!(■  - J.f)  is  the  unit  step 
function.  With  this,  the  flow  tanqency  condition  becomes 

^ = f‘(.:.)  - a (r)/6  - H(f.  - ?if)fl(T)/<S  on  n 0 for  0 i !!  i 1.  (12) 

We  note  that  this  condition  is  applied  on  n * 0 as  is  consistent 
with  the  small  disturbance  assumption.  In  addition,  to  the  ordei 
considered  in  equation  (2°*)  there  is  no  explicit  dependence  on  a 
in  the  surface  boundary  condition  (12).  The  wake  condition  is 
applied  as 

[f . ] = 0 on  n = 0 for  ,f,  x 1 ( i 1 ) 

where  the  square  brackets  indicate  the  "jump"  in  the  enclosed 
quantity  and  the  condition  is  again  applied  alono  n - 0.  The 
far  field  boundary  condition  is 

Mu)"  + (if^)  0 as  f*  + n‘  -a».  (3A) 

so  that  equation  (29)  and  correspond i nq  boundary  conditions 
(12) -(U)  complete  a description  of  the  flow  field  problem  if  an 
initial  profile,  <M-r.,n,o),  is  specified. 

With  this  formulation,  the  unsteady  pressure  coefficient  is 
given  by 


2/J 

C = -2<S  * . 

P V,  * 


(15) 


The  correspond i nq  lift  and  moment  coefficients  follow  as 


c»  - lo  V * Vdr- 


(If) 
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SECTION  ITT 

THE  METHOD  OF  SOLUTION 


We  now  seek  a method  of  obtaininci  solutions  to  the  coupled 
system,  equation  (261  with  appropriate  initial  conditions  and 
equation  (2*))  with  boundary  conditions  (.12) -(14)  and  initial 
profile  , n,o)  . This  formulation  wall  be  considered  qeneral 
with  respect  to  equation  (26)  such  that  anv  number  of  decrees  of 
freedom  are  permitted.  Equation  (261  is  first  written  as  an 
equivalent  first  order  system  by  definino 

x , ( r ) = x ( t ) 


and 


X • ( 1 1 = X ' ( t 1 


(•’'11 


so  that 


X { ( T ) X : ( t 1 


(•Ml 


and  , 

X:'(r)  = M~  [C(i)/y  - PX,  (t)  - A.X,(t)].  (•'') 

Fouations  (41)  and  (421  were  then  solved  bv  a 4 -point  Adams-Moul t on 
inteqration  in  time  (201.  This  is  an  implicit  pred i ctor-cor rector 
technioue  havino  a local  truncation  error  of  order  (a  t l . a 
number  of  other  methods,  both  explicit  and  implicit,  were  tried. 

It  was  found  that  when  coupled  with  the  aerodynamic  equation, 
calculations  eventual lv  became  unstable  for  central  difference 
and  i-term  Taylor  series  explicit  schemes  and  for  Fulor-i’auohv 
and  Milne  1-point  implicit  schemes,  even  when  extremely  small  time 
steps  were  employed.  A Milne  5-point  implicit  technioue  with  local 

l l 


truncation  error  of  order  (Ax  ) , although  stable,  appeared  to 

produce  less  accuracy  than  the  Adams-Moulton  integration. 

It  will  be  of  interest  to  compare  some  time  integrated  results 
with  homogeneous  solutions  to  the  structural  system.  In  particular, 
for  , the  homocieneous  result  should  be  recovered.  For  the 
case  of  C(t)  = 0 the  associated  eigenvalue  problem  correspond  inn 
to  equations  (41)  and  (42)  was  considered  and  the  eigenvalues 
and  eigenvectors  obtained.  The  solution  vector  was  then  represented 
as  the  linear  combination  of  the  eigenvectors  which  satisfied  the 
specified  initial  conditions.  These  solutions  are  thus  considered 
"exact"  as  they  are  free  of  truncation  error  which  is  present  in 
time  integrated  results.  No  attempt  was  made  to  obtain  normal 
modes  as  this  is  not  in  general  possible  for  arbitrary  values  of 
the  mass,  damping,  and  stiffness  matrices  (M,P,K)(21). 

An  efficient  time-implicit  finite  difference  algorithm  (13) 
has  been  developed  for  obtaining  solutions  to  the  low  frequency 
transonic  equation.  This  technique  was  incorporated  in  the  computer 
code  LTRAN2 , developed  by  Ballhaus  and  Oooriian  (16)  for  the 
purpose  of  computing  unsteady  transonic  flows  over  airfoils 
using  the  form  of  the  potential  equation  given  by  (29)  . Several 
minor  changes  in  the  basic  code  were  necessary  in  order  to 
accommodate  a simultaneous  solution  of  the  structural  equations 
(41)  and  (42).  Details  of  LTRAN2  are  briefly  summarized  here. 


The  basic  LTRAN2  code  employs  a non-iterative  alternating- 
direction  implicit  (ADI)  scheme  to  advance  the  solution  for  the 
perturbation  potential,  4* , from  one  time  level  to  the  next  at 
each  grid  point  in  the  computational  flow  field.  Differencing 
in  the  :.  -direction  is  of  the  mixed  type  which  has  been  quite 
successful  in  maintaining  stability  for  both  subsonic  and 
supersonic  flow  regions.  Conservation  form  of  the  equation  is 
preserved,  which  is  essential  for  a proper  description  of  shock 
wave  motions.  While  the  ADI  scheme  has  no  time  step  limitation 
for  stability  based  upon  classical  linear  stability  analysis, 
instabilities  may  be  generated  by  the  motion  of  shock  waves  due 
to  the  mixed  differencing.  Thus,  At  must  be  chosen  such  that  shock 
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waves  do  not  travel  more  than  one  mesh  point  in  the  : - direction 
over  a single  time  step.  A cubic  spline  is  used  to  approximate 
the  airfoil  geometric  function  f(.\),  which  is  usually  provided 
as  tabular  data  for  NACA  airfoils.  Aerodynamic  moments  are 
evaluated  by  Simpson's  rule  integration  according  to  equations 
(if)  - (IS).  A smooth  non-uniform  computat ional  mesh  which  is 
symmetric  about  n = 0 is  employed.  The  grid  spacing  is  such  that 
points  are  clustered  near  the  airfoil  leading  and  trailing  edges 
in  the  i,  - direction,  and  near  n = 0 in  the  n - direction.  Mesh 
boundaries  are  taken  sufficiently  far  from  the  airfoil  such  that 
condition  (14)  may  be  approximated.  The  airfoil  surface  is 
described  by  if  mesh  points.  Details  of  the  grid  system  are  the 
foil ow i ng : 


a . 

number  of 

- iv<  i n t s 0 o 

b. 

number  of 

n-oo ints  7R 

O • 

\ •.  . 0 . 

min 

00  3 30 

d. 

An  = o . 

min 

02000 

e . 

-1033.5304 

7 x t,  ^ 855.0111 

f 

-811. 12200 

» • 811.12200 

Results  of  computations  from  I.TRAN2  have  been  shown  to 
compare  well  with  solutions  of  the  time-dependent  Kuler  equations, 
an.)  have  reproduced  unsteady  transonic  behavior  which  has 
commonly  been  observed  experimentally  ih'' . In  ordet  to  accommodate 
the  simultaneous  solution  of  the  structural  equations  (411  and 
(4.')  it  was  necessary  to  modify  l, TRAN.'  such  that  the  time  marching 
technique  is  now  a two-step  (pred i ctor-correct oil  procedure. 

Although  this  in  effect  doubles  the  computing  time  for  a tixed 
\ : , without  this  modification  stable  calculations  could  not  in 
general  be  obtained. 


SECTION  IV 
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RESULTS 

In  this  section  the  behavior  of  numerical  solutions  obtained 
by  th^>  method  outlined  in  Section  7.0  is  considered.  Accuracy 
of  the  solutions  is  discussed  and  the  choice  of  structural 
parameters  indicated.  Results  are  presented  for  an  airfoil 
having  both  one  and  three  deqrees  of  freedom.  Solutions  were 
obtained  for  fixed  initial  conditions  and  several  values  of  the 
reduced  density, y . Additional  results  indicate  that  for  y fixed, 
stable  or  unstable  motion  can  be  produced  solelv  by  the  choice 
of  initial  conditions.  All  of  the  calculations  presented  here 
correspond  to  an  NACA  64A010  airfoil. 

4.1  Details  of  the  Computations  and  Accuracy  of  the  Solutions 
In  order  to  indicate  the  stability  of  the  computational 
method,  it  will  be  demonstrated  that  an  exact  solution  of  the 
coupled  aeroelastic  system  of  eauations  can  be  reliably  reproduced 
by  time  integration.  For  this  purpose,  the  airfoil  motion  was 
forced  for  three  periods  of  oscillation  according  to  the  prescribed 
functions . 

o(r)  = ct  ( t ) = B ( t ) = 0.01745  sin  ( 4 . 3t  ) . (43) 

Using  a steady-state  profile  as  the  initial  condition,  equation 
(29)  was  inteqrated  in  time  bv  I, TRAN? . After  a short  period  of 
time  the  effect  of  the  initial  conditions  became  nealigible  such 
that  the  aerodynamic  forces,  C\,C  , and  C . were  periodic. 
Analytical  expressions  for  these  forces  were  then  extracted  from 
numerical  results  and  introduced  into  equations  (41)  and  (42). 

The  structural  parameters  comprisina  the  mass,  damoinq,  and  stiff- 
ness matrices  (M,P,K)  and  the  reduced  density  y were  then  chosen 
such  that  the  analytic  solution  to  equations  (41)  and  (42)  was 
aiven  by  the  previously  prescribed  functions  (4  7)  for  forced 
motion.  At  this  point  the  coupled  system  was  inteqrated  in  time 
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using  as  initial  values  the  conditions  present  when  the  forced 
motion  was  terminated. 

This  procedure  was  used  to  study  both  stability  and  accuracy 
characteristics  of  various  numerical  techninues  for  intearatinn 
the  structural  equations  of  motion.  It  was  found  that  Cf  was  the 
aerodynamic  coefficient  most  sensitive  to  procedural  variations. 
One  result  of  this  studv  is  shown  in  Finure  2 for  the  case 
M = 0.72.  The  time  dependent  lift  is  displayed  for  both  forced 
and  free  motion.  Over  three  periods  of  free  oscillation  tne 
variation  in  c,  was  never  more  than  0.5%  and  it  appeared  that  the 
inteoration  could  have  proceeded  indefinitely  without  instability 
and  with  reasonable  accuracy. 

We  now  devote  our  attention  to  solutions  of  the  coupled 
aeroelastic  equations  for  representative  choices  of  the  structural 
parameters.  Impulsive  motion  from  an  equilibrium  state  will  be 
the  only  tvpe  of  motion  considered  here.  For  all  the  results 
presented  in  this  section,  initial  values  correspond  to 

a (0)  = a (0)  = 6(0)  = 0, 

o ' (0)  = P ' (0)  = 0, 
a ' (0 ) prescribed 

with  <ME,n,o)  qiven  bv  the  steadv-state  profile  for  Mot  = 0.82. 
The  initial  pressure  distribution  for  these  results  is  shown  in 
Ficure  2.  It  is  noted  that  this  case  is  supercritical  with  a 
shock  appearinq  between  r,  = 0.50  and  E = 0.60.  It  should  be 
pointed  out  that  solutions  are  not  limited  to  the  above  initial 
values,  however  it  was  necessary  to  fix  a certain  number  of 
parameters  in  order  to  conduct  a reasonable  study.  Due  to 
inherent  nonlinearity  in  the  transonic  equation,  various  types 
of  motion  may  result  dependinq  on  the  choice  of  a ' (o)  . This  is 
quite  unlike  classical  flutter  analysis.  Physically,  the  initial 
conditions  o'(o),  a ' (o) , and  P ' (o)  may  be  interpreted  as  impulsive 
qusts  striking  the  airfoil.  By  varvinq  these  initial  conditions, 
stability  boundaries  may  be  established.  Combinations  of  these 
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FIGURE  2.  Comparison  of  Forced  and  Free  Oscillations 
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conditions  need  necessarily  be  considered,  however,  as  super- 
position cannot  be  applied. 

The  airfoil  with  pitch  point  located  at  the  quarter  chord 
aileron  leading  edge  at  75%  of  the  chord,  and  aileron  hinge  point 
at  80%  of  the  chord  is  depicted  in  Figure  1.  These  choices 
correspond  to  = 0.25,  = 0.75  and  = 0.80.  Structural 

parameters  were  selected  as  the  following: 

?,  = -0.1818, 
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Various  values  of  the  damping  coefficients  ( r. 


and 


, 1 and 


a or 

the  reduced  density  ( u ) were  chosen  and  will  be  noted  in  specific 
results.  By  way  of  comnarison,  if  the  airfoil  is  assumed  to  consist 
of  a homogenous  material  of  uniform  densitv,  then  the  following 
mav  be  obtained: 

l0  = -0.1818, 


r - 

'•f 


and 


f*  = 

'o 

. * 

‘•f 


0.0034  , 
0.0846, 
= 0.0052. 


For  all  calculations  the  time  step,  At,  was  specified  to  be 
0.01745.  This  corresponds  to  360  time  steps  per  period  of 
oscillation  at  a reduced  freguencv  (u>c/Uv)  of  0.215.  Time 
accuracy  of  the  solutions  was  further  confirmed  by  doubling  this 
nominal  value  of  At  and  comparing  with  previous  results.  Tt  was 
found  that  the  solutions  agreed  to  three  significant  figures  in 
all  dependent  variables,  over  several  thousand  time  steps  even  for 
cases  which  exhibited  unstable  motion.  When  the  time  increment 
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was  increased  by  a factor  of  five,  however,  the  computations 
eventually  became  unstable.  This  was  most  likelv  due  to  a 
violation  of  the  time  step  size  limitation  with  respect  to  the 
motion  of  the  shock  across  the  computational  mesh  which  was 
previously  described.  All  calculations  were  performed  on  a 
CDC  Cvber  74  computer  and  required  approximately  20  minutes  of 
computinq  time  for  every  1000  time  steps. 

4.2  The  One  Decree  of  Freedom  Case 

It  is  useful  to  study  the  behavior  of  solutions  for  a 
sinqle  degree  of  freedom  airfoil.  Recause  there  are  no  super- 
imposed modes  of  motion  in  this  case,  the  characteristic  details 
of  the  solution  are  quite  clear.  For  an  airfoil  oscillatinq  in 
pitch  only,  the  structural  system  reduces  to  equation  (24)  with 
0 = 6=0.  Fiqure  4 indicates  the  time  resnonse  of  the  sinale 

degree  of  freedom  airfoil  for  £ = 0.01  and  a'(o)  = 1 . 0 for 
two  values  of  the  reduced  density.  The  homogeneous  solution  is 
also  shown  for  comparison.  Here  the  pitching  displacement  has 
been  normalized  by  the  impulsive  pitch  velocity  such  that  the 
homogeneous  solution  is  universal  for  all  choices  of  a'(o)  . Thi 

will  prove  a convenient  form  to  use  for  later  results  when  varia- 
tions m the  initial  condition  are  considered.  The  time  axis  is 
represented  in  terms  of  the  number  of  time  steps,  N,  in  order  to 
emphasize  the  fact  that  the  resolution  of  the  discretization  in 
time  is  sufficient  to  consider  the  time  history  a continuum.  The 
corresponding  physical  time  is  easily  obtained  from  equation  (10) 

t = cAtiS  N/u  . (44) 

00 

For  the  choices  of  At  and  5 considered  here  this  results  in  one 
chord  length  of  airfoil  travel  for  every  12.344  time  steps. 

With  = 0.01  it  is  seen  that  the  homogeneous  solution  is 

very  slightly  damped.  For  u = 50,  the  time  integrated  solution 
varies  slightly  from  the  homogeneous  result.  In  fact  for  ^ = 100 
the  homogeneous  solution  was  recovered  virtually  intact.  Once 
again  this  confirms  accuracy  of  the  numerical  method.  As  the 


reduced  is  increased,  a shift  in  the  oscillation  frequency  as  well 
as  well  as  an  increase  in  amplitude  is  noted.  Corresponding  time 
histories  for  the  moment  coefficient  about  the  pitch  point  are 
displayed  in  Figure  S,  and  are  seen  to  exhibit  a behavior  similar 
to  that  of  the  pitching  displacement. 

For  u - 100  and  • 0.01  the  pitchinq  displacement  for 

a 

several  values  of  the  initial  condition  is  shown  in  Figure  f> . 

The  case  ot  (o)  = 4.0  is  nracticallv  neutrallv  stable.  Tt  is 
seen  that  for  a ' (o)  = 1.0  the  amplitude  of  oscillation  is  about 

JO*  damped  in  two  periods  of  oscillation.  Slightly  growing 
amplitudes  are  produced  with  a * (o)  = 7.5.  This  behavior  exempl  i- 

fies the  nonlinear  character  of  the  coupled  aeroelastic  system. 

Tt  is  noted  that  no  appreciable  change  in  frequency  is  evident 
for  variation  in  the  initial  conditions.  Moment  coefficients 
corresponding  to  these  cases  appear  in  Figure  7.  For  small  time 
(N  < 200),  the  anomalous  behavior  apparent  in  the  cases  d (o)  =4.0 
and  n (o)  = 7.5  is  most  likely  due  to  starting  phenomena  and  an 
adjustment  in  the  phase  difference  between  a and  C 

mo 

4.3  The  Three  Degree  of  Freedom  Case 

We  now  focus  our  attention  on  the  three  degree  of  free- 
dom airfoil.  For  the  results  considered  in  this  section  all 
damping  coefficients  were  assigned  the  identical  value  of  0.01 
(i.e.,  t = Cp  = 0.03).  Unsteady  displacements  for  the 

initial  condition  n‘(0)  = 1 and  several  values  of  reduced 
density  are  given  in  Figures  8-10  with  the  homogeneous  solution 
indicated  for  comparison.  The  plunging  displacement  shown  in 
Figure  8 is  seen  to  vary  from  the  homogeneous  result  and  experience 
a shift  in  frequency  as  u increases.  Damped  motion  appears  to 
occur  for  u = 100  and  p = 50,  but  an  increasing  amplitude  is 
noted  for  u = 25.  In  Figure  8 this  same  behavior  is  observed  for 
the  pitching  displacement.  There  seems  to  be  very  little  snift 
in  frequency  for  the  cases  p = 100  and  p = 50.  It  appears  that 
for  large  time  the  response  is  tending  to  make  the  plunge  and 
pitch  frequencies  identical.  Figure  10  indicates  the  aileron 
displacement  for  these  cases.  It  is  seen  that  there  is  apparent Iv 
damped  motion  with  very  little  frequency  shift  for  all  values  of  . 
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Corresponding  aerodynamic  coefficients  for  these  cases  are 
found  in  Figure  11-13.  In  general,  they  exhibit  behavior  which 
is  very  similar  to  that  of  the  displacements.  The  lift  and 
pitching  moment  are  damped  for  u = 100  and  u = 50,  but  are 
growing  with  U = 25.  For  ail  values  of  u the  aileron  moment 
is  damped.  It  is  again  indicated  that  for  u = 25  the  motion 
consists  primarily  of  a coupling  between  the  plunge  and  pitch 
modes  onlv. 


As  in  the  case  of  the  one  degree  of  freedom  airfoil,  we 

now  consider  the  time  response  for  fixed  u and  several  choices 

of  the  initial  condition,  a*(0).  Unsteady  displacements  for 

u = 25  appear  in  Figures  14-16.  Time  integrated  results  are 

now  compared  with  the  undamped  homogeneous  solution 

(i.e.,  r = t = C 0 =0).  It  has  alreadv  been  indicated  that  with 
oat 

U = 25  growing  amplitudes  of  the  displ acements  occurred  when 
compared  with  the  damped  homogeneous  result.  But  as  the  damped 
homogeneous  solution  is  itself  decaying  this  may  be  mi s lead i no, 
especial lv  because  of  the  several  superimposed  modes  which  are 
present  m the  solution  of  any  one  of  the  displacements.  The 
long  time  behavior  can  alwavs  be  deduced  by  integrating  suffi- 
ciently far  in  time.  If  we  wish  onlv  to  establish  whether 
oscillations  are  crowing  or  decaying  bv  observino  the  relativelv 
small  time  response,  then  comparison  with  the  undamped  homogeneous 
result  may  prove  quite  useful . 

The  unsteady  pitching  displacement  is  shown  in  Figure  14. 

While  ail  solutions  have  shifted  in  frequencv  from  the  homogeneous 
res.ult,  little  variation  in  frequency  occurs  for  chances  in  the 
initial  condition.  All  results  appear  to  have  growing  amplitudes, 
but  less  so  for  a ( 0 ) = 5 then  for  a’(0)  = 1 and  n‘(0)  = 3. 

The  trend  with  respect  to  frequency  is  verv  similar  for  the  pitching 
displacement  in  Figure  15.  Here,  however,  the  solution  for 
a’{0)  = 5 is  clearly  growing  whereas  the  other  results  mav  in 
fact  be  damped.  This  is  also  true  for  the  aileron  displacement 
shown  in  Figure  16. 
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Fiqures  17-10  display  the  correspondinq  solutions  for  the 
aerodynamic  coefficients  with  u = 25.  The  lift  and  pitchinn 
moment  behavior  is  quite  like  that  of  the  displacements.  Little 
variation  in  frequency  is  observed  for  chances  in  the  initial 
condition  and  qrowinq  amplitudes  are  clearly  indicated  for  the 
case  a (0)  = 5.  The  abrupt  chancres  in  the  slope  of  the  aileron 
moment  are  due  to  the  unsteady  shock  wave  oscil latino  across  the 
control  surface.  This  is  evident  from  the  surface  pressure 
distribution  shown  in  Fiqure  20  for  N = 115.  The  time  here 
corresponds  to  relative  maximums  in  the  itchinn  and  aileron 
moments  and  a relative  minimum  in  the  lift.  For  u'(0)=  5, 
the  entire  upper  surface  is  now  subcritical  while  most  of  the 
lower  surface  has  become  supercritical  with  the  shock  located 
near  the  trailinq  edqe.  Recall  that  in  the  initial  profile  a 
shock  was  present  on  both  surfaces  near  mid-chord  (see  Fiqure  O. 
Ry  comparison,  for  the  case  «'(0)  = 1 the  shock  is  displaced 
only  sliqhtly  forward  on  the  upper  surface  and  sliqhtlv  aft  on 
the  lower  surface  from  its  initial  position.  This  type  of  shock 
wave  motion  is  representative  of  that  which  can  be  simulated  bv 


the  LTRAN2  code. 


To  establish  the  lonq-time  behavior  of  the  solutions,  the 
case  a’(0)  = 5 and  u = 25  was  inteqrated  in  time  for  several 


thousand  time  steps.  The  extended  time  histories  of  the  dis- 
placements for  this  computation  are  shown  in  Fiqure  21.  As  was 
predicted  from  the  short-time  behavior,  all  displacements  clearlv 
exhibit  increasinq  amplitudes.  The  pitchinq  and  aileron  dis- 
placements are  oscillatinq  at  approximately  the  same  frequency 
which  is  about  twice  that  of  the  plunqinq  displacement.  For  the 
aerodynamic  coefficients  shown  in  Fiqure  22,  however,  all  appear 
to  have  the  same  oscillation  frequency.  This  extended  time  history 
indicates  that  such  calculations  may  be  carried  out  for  marqinally 


stable  cases  when  the  behavior  cannot  be  deduced  from  a small-time 


solution.  The  lonq-time  solution  may  in  qeneral  be  questionable 
due  to  the  accumulation  of  truncation  error  in  the  time  inteqration, 
Results  of  the  forced  oscillation  study,  however,  indicate  that  the 
solutions  presented  here  can  be  considered  reliable. 
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Figure  21.  Extended  Time  History  of  Unsteady  Displacements  for  Three  Degree 
Freedom  Airfoil  with  * ' (0)  = 5 and  u = 25. 


Figure  22.  Extended  Time  History  of  Unsteady  Aerodynamic  Coefficients  for 
Three  Degree  of  Freedom  Airfoil  with  a'(0)  = 5 and  u = 25. 


SFCTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

A procedure  has  been  developed  for  the  aeroelastic  analvsis 
of  a two-dimensional  airfoil  in  transonic  flow.  The  method  has 
been  shown  to  be  stable  and  fv^curate.  Solutions  have  been 
presented  for  both  one  and  three  decree  of  freedom  airfoils 
usinc  representative  values  of  structural  parameters.  Stable 
or  unstable  modes  of  oscillation  were  shown  to  depend  upon  the 
choice  of  initial  values. 

An  obvious  improvement  of  the  work  presented  here  is  the 
extension  to  three  spatial  dimensions.  At  the  present  time,  the 
onlv  efficient  method  of  solvinc  the  three  dimensional  fluid 
dynamic  equations  for  the  unsteady  transonic  case  is  bv  harmonic 
analvsis  (3,  5,  & 6)  which  eliminates  the  nonlinear  time  dependent 
nature  of  the  flow.  It  is  felt  that  in  the  near  future,  however, 
time  inteqration  of  the'  small  disturbance  potential  equation  in 
three-dimensions  mav  be  accomplished  with  reasonable  oomputina 
times.  The  airfoil  can  then  be  represented  as  a flexible  cantilever 
beam  by  a direct  extension  of  the  procedure  outlined  in  this  report. 
Such  an  analysts  would  clearly  provide  a treatment  which  more 
closely  resembles  the  physical  situation. 

Due  to  the  assumption  of  low  frequency  oscil  lat ions , the  term 
••  does  not  appear  in  the  form  of  the  potential  equation  (2°)  used 
in  this  report.  While  it  is  felt  that  in  neneral  this  assumption 
vs  iustified,  some  question  must  be  raised  about  the  validity  of 
the  very  small-time  behavior  of  solutions,  particularly  for  the 
cases  of  impulsive  motion  considered  here.  Nealect  of  this  term 
probably  has  little  effect  on  the  lona-time  airfoil  response.  At 
the  present  time,  no  efficient  fullv  time  implicit  scheme  exists 
for  inteqration  of  the  potential  equation  witti  the  term  included. 

Finallv,  the  effects  of  viscosity  should  not  be  neglected 
if  aerodynamic  coefficients  are  to  be  accurately  predicted. 
Instantaneous  shock  locations,  and  hence  resultant  moments,  can 


42 


u 


depend  greatly  upon  unsteady  viscous  interaction.  Couolinq  of 
the  viscous  effects  with  inviscid  flow  field  calculations  is  not 
a simple  task.  In  the  case  of  unsteady  transonic  flows  only  one 
procedure  for  accomplishina  this  has  been  suaqested  (9) , and 
its  present  application  remains  largely  an  art. 
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